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<^ ; Abstract 

^ ■ In this letter a general method for expanding paraxial beams into multipolar electromagnetic 

t^ . fields is presented. This method is applied to the expansion of paraxial modes with orbital angular 
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momentum (0AM), showing how the paraxial 0AM is related to the general angular momentum of 
an electromagnetic wave. This method can be extended to quasi-paraxial beams, i.e. highly focused 
laser beams. Some applications to the control of electronic transitions in atoms are discussed. 



The angular momentum of light fields continues to draw the interest of the scientific 
community. Since the work of Allen and coworkers, showing that the angular momentum 
components of a paraxial light beam, i.e. the spin and the orbital angular momentum 
(0AM), could be controlled separately [l|], there has been an ongoing discussion regarding 
concepts as the conservation of angular momentum in different light interaction processes 

0.0. 

hght 



4|, the angular momentum flux |5| or the measurement of the angular momentum of 

sa. 



In parallel to this scientiflc discussion, there has been an explosion of applications of the 
angular momentum of light, and in particular to the orbital part, covering areas as diverse as 



111 ] or quantum 



trapping and rotating microparticles 8|, astrophysical measurements 

information 121. Also, recently there has been a renewed interest in the interaction of the 

n 

orbital angular momentum of light with atomic ensembles IIJ and the electronic structure 

n 

of atoms and molecules [IJ] . 

On the other hand, a general electromagnetic fleld it is known to have an angular mo- 
mentum content 15| . One can even construct solutions to the electromagnetic fleld with the 
right spherical symmetries to have well deflned values of the total angular momentum and 
one component of the angular momentum. This kind of solutions receive the general name 
of multipolar solutions of the electromagnetic fleld. Up to now, the multipolar flelds and 
the paraxial modes with 0AM have been used separately in different regimes. In this letter 
I try to close this gap and I calculate the multipolar content of any paraxial beam and in 
particular of the 0AM modes. 

A general electromagnetic fleld in free space has two components of the angular momen- 
tum (J): the spin part (hereby denoted by S), which is related to the vectorial character 
of the fleld, and the orbital angular momentum (0AM) L, which is related to the spa- 



tial structure of the fleld 



15| . Both parts are needed to generate rotations in space, and 



consequently only the combination of the two components plays a meaningful role in the 
rotational symmetries of an electromagnetic wave, i.e. J = L + S [16|. Multipolar modes are 
precisely a set of solutions to the Maxwell equations which are eigenvectors of the square of 
the total angular momentum J^, one component of the total angular momentum, i.e. the 
^-component, J^ = L^ + S^ and the parity operator. The exact form of these flelds can be 



found for example in 16|. In this letter, I will just use the properties of the vector potential 
of the multipolar monocromatic flelds in the solenoidal gauge. Then following ref. [la] the 



vector potentials of the set of multipolar fields is A^-^, where x = m represents the magnetic 
multipoles and x = e the electric multipoles. Both classes of multipoles have eigenvalues 
3'^AjJ^ = j{j + l)Aj-^ and Jz-^jm = "^A^^, but magnetic and electric multipoles differ in 
their parity. Any given electromagnetic field can be fully expanded onto multipolar fields, 
if we also add a longitudinal field of the form Ajl^ = V^jYjm, with Yjm being the usual 
spherical harmonics. 

Let me start by constructing an electromagnetic field by superposing rotated circularly 
polarized plane waves: 

/■TT f'2'K 

A= / sm{9)d9 d^g(9,^)R{6,(p)epexp{ikz) (1) 
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where the original plane wave propagates along the z direction and is right (left) handed 
polarized when p = 1 {p = —1). The altitude and azimuthal rotation angles are given 
respectively by {O^ip). The operator i?, rotates the vector field in the usual manner: 
R{6,ip)epexp{ikz) = M~^epexp{ikz ■ (Mr)). Note that this field is completely general 
and does not have to fulfil paraxiality. In particular, it can be used to describe beams highly 
focused with aplanatic lenses [l7|. The function g{9,(f) controls the relative amplitudes of 
the plane waves. 

The main idea of this letter is to find a family of fields using suitable functions g{9, if), 
which in the limit of paraxiality can be identified with the usual orbital angular momentum 
modes. Let me then develop the rotation operator to the lowest orders in the angle 9 to 
obtain 
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where I use cylindrical coordinates (p, 0,2;). Now, compare this result with a circu- 
larly polarized paraxial field of the kind A = Fi{p, z) exp{ikzz) exp{il(f))ep. Any parax- 



ial beam can be expressed as a convenient superposition of this kind of modes \iB\, so 
our results apply to any paraxial beam. Besides, this kind of beams are cylindrically 
symmetric and thus are eigenmodes of both the 0AM and the spin AM components in 
the z direction, i.e. L^ and Sz- Note, however, that this feature is in principle only 
valid in the paraxial approximation, when kz ^ k and Fi is slowly varying in z with 
respect to 1/k. This paraxial field can be reexpressed in its Fourier components as 
A = /o°° kr dkr /o'^ d(f)kfi{kr) exp(i/0fc) exp{ikrrcos{(f) — 0fc)) exp{ikz — ikl/{2k)z)ep. 



It is clear from inspection of Eqs. ([T]) and ([2]) that if we do the correspondence 

kr = ksm{9), (pk = V^, 

g{kr/k, (f)k) = k"^ fi{kr) exp {i{l + p)0fc) , (3) 

then the two expansions of the paraxial beam, i.e. in Fourier components and in rotated 
plane waves Eq. ([T]), are one and the same, within the paraxial approximation. Then, 
the set of functions g, fulfilling Eq. ([3] defines the sought-after family of fields. Note that, 
in general, we can use functions fi{kr) which are not paraxial, thus defining cylindrically 
symmetric electromagnetic fields. 

This identification is useful, because now 1 will use the expansion of a rotated plane wave 
into multipole waves [l6|, [l9| and perform the integral over the rotation parameters. The 
expansion reads 

oo j 

Di^{^,eM^'^;^ + wA-^^l (4) 

where Di^p(if, 6, 0) are the matrices of rotation for irreducible tensors of order j. This matrix 
can be expressed as D-i^p{(p,9,0) = exp{—ii7np)d{^p{9), where 

dLM = lU + p)\{j - p)\{j + m)\{j - m)\f' X 
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cos 7^1 (-sin-) • (5) 



If we put together Eqs. ([T]),(ll]) and ([5]), we obtain 

oo j „7r 

A = E E i'{'^j + if" j MO)dedipie) 

j=l m=-j 

J dipexp{-imip)g{9,ip) [aJ™^ + ^aJ™] . (6) 

Now, I will use the cylindrically symmetric function g from Eq. (^ to perform the 
integrals. In this way, the ip integral is immediate and the following result holds 

oo 

A = 2vr E ^C^J + ^y^'C,JA^l.+^A^^ 



j=\\i+p\\ 



■^i (i+p) "^ "^j {i+p). 



c,i, = J^ sm{e)dedi,^p^p{e)Mkszn{e)). (7) 



This equation is the main result of this letter. Note that Eq. ([7]) is a valid solution of the 
Maxwell equations, as fi{kr) is not restricted to paraxial beams, and then can be applied to 
a wide number of situations. Also, note that all the multipolar beams in the superposition 
share the same value oi m = I -\- p. This has the obvious meaning of summing the 0AM 
and SAM in the paraxial approximation, but in the more general case, simply implies that 
we have been able to form a family of fields with a well defined Jz value. This solution is 
different from that presented in Refs. l|, |5|, l20(] and it does not have to fulfil the angular 
momentum equations derived there. 

On the other hand, this result has an interesting side effect. As multipole fields only 
depend on the wavelength, we can control the multipolar content of a field with only paraxial 
fields. I will step back a little bit and assume the paraxial approximation again. In this 
case, the amplitude of the multipolar field of order (j, m = I + p) can be approximated by 

Cjip = {-y[u+p)Kj-py.{j + i+p)Kj-i-py.f' 

^ {j - I - P - s)\{j + p - s)\{s + ly.sl 2'+2« ^^ 

where the number Ma{fi} = J^ K dky.k'^ fi{kr) is the momentum of order a of the function 

/^ 

Let me now discuss an important class of functions, the Laguerre-Gaussian modes (LG). 
The LG modes are defined by two indices, the 0AM index, which I conveniently call / and 
an index stating the number of transversal nodes of the function, which I call q. Then, the 
fi functions read fi^,{K) = ^^ ^^^{nJoK I V^f^ L{l(w;2fc2/2) e^wiejA + i{q - |/|/2)7r), 
where wq is the beam width of the LG mode in real space, and L^q{x) are Laguerre polynomi- 
als. When this expression is used, the set of momenta Ma can be calculated analytically and 
for the special case of LG modes with g = 0, we can write the sum in Eq.® in a strikingly 
compact form 
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{j+py-{j-p)Kj+i+py-{j-i-py- 
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I will detail now some examples. The simplest case is that of a paraxial Gaussian beam 
{wq ^ 1/A;) with circular polarization. In this case, formula flTOj) gives the following result 

oo 

A = 2vr 5: z^-(2j + ly/'C, ,, [aJ™) + ipA^ 
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FIG. 1: Multipolar decomposition of a Gaussian beam. First colunm: Right circular polarization 
(p = +1). Second column: Left circular polarization (p = —1). First row: weight of the multipolar 
components C?q . Second, third and last rows represent the different components (respectively 
right circular, longitudinal and left circular) of the field at 2 = 0, after performing the summation 
of the multipole components. Note the different scales for the different components. 
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Note the following features of this field. First, the field is composed of multipolar solutions 
with different total angular momentum, but with the same projection of the angular mo- 
mentum in the direction z, which is +1 or —1 depending on the polarization of the field. 
Also, in both cases the amplitudes of the multipolar components are the same, as Cj-op does 
not depend on the sign of p. In Fig. ([1]) we plot the different amplitudes of the multipolar 
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FIG. 2: Multipolar decomposition of Laguerre Gaussian beams. First column: Right circular 
polarization {p = +1) and azimuthal index negative (/ = —1). Second column: Left circular 
polarization (p = —1) and azimuthal index positive (/ = +1). First row: weight of the multipolar 
components: C'^±i_i- In this case Cj +i-i = —Cj _i+i. Second, third and last rows as in Fig. 1. 



components for Wq = 15/ k, and the field resulting when sum fllOl) is numerically performed. 
Another interesting case occurs when we use Laguerre-Gaussian beam with / = ±1, g = 
and circularly polarized p = —I. I have shown that, in the general case, all the multipolar 
fields have a component of J^ is equal to m = I + p. Then, in this example the projection of 
the angular momentum in the z direction of both considered fields is zero, i.e. m = l+p = 0. 
The multipolar decomposition gives the following result. 
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2V2+i^-2^L]_, {{wo/kr) (11) 

Note that the weights of the multipoles of the two fields have the same magnitude and 
opposite signs: Cj +i_i = —Cj +i_i. An example of such fields is found in Fig. ([2]). 

This particular example is very interesting as the two considered fields are decomposed in 
the same subset of multipolar modes: Aj'^ and Ajg. This will always be the case of fields 
with opposite polarizations and A/ = 2. However, in this particular example, both fields 
share the same weights in the decomposition. This means that if we produce a superposition 
of the kind: A = aLGiQe^i + jSLGioe^i, we can produce a purely transversal electric field 
when a = —j3 or a pure magnetic one with a = f3. In other cases with opposite polarization 
and A/ = 2, one can produce a multipole field with a well defined parity just in a set of 
discrete values of j, by playing with the amplitudes of the superpositions and the beam 
widths of the two fields. 

The examples above, and the formulae ([7]) and (TTOl) demonstrate the control of the su- 
perpositions of a set of multipolar fields by using paraxial beams. As considered above, 
multipolar fields are the solutions of the electromagnetic field one should use when treating 
spherically symmetric problems. For example, those problems where there is an exchange of 
angular momentum between the electromagnetic field and material particles are particularly 
well suited to be treated with multipolar expansions of the electromagnetic field. One such 
case are the electronic transitions in atoms or molecules. It is well known that multipolar 
transitions are increasingly more difficult to excite (or de-excite) for larger j's. This is why 
almost all the literature dedicated to the exchange of angular momentum (or more exactly 
orbital angular momentum), treats the problem in the dipolar approximation. Actually an 
easy calculation shows that the probability for a multipolar transition to occur scales as 
{a/Xy^^ , with a being the typical radius of the system involved, and Aj the multipolar 



transition 151 ]. 



Here we are interested in a different, but related problem. Let's consider the case where 
it is needed to control a certain multipolar transition. This is the case, for example, when 
our system is in certain metastable states, as is the case of trapped Ca"*" ions in quantum 
information applications (21|. In those cases the rate of transition is low compared to dipole 
transitions, simply due to the small size of the ion, compared with the wavelength of the 
field. Nevertheless, we can still maximize the overlap of the laser beam we are using to induce 
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the transition with the electromagnetic multipole field associated with this transition. Also, 
in practical applications we have to deal with paraxial or near-paraxial beams. In this 
situations we can use Eq. (jS]) to engineer the shape and polarization of our control beam, 
to maximize the overlap with the desired transition. 

Importantly, our method allows to control not only the total angular momentum exchange 
involved in a certain transition, but also the exchange of Jz- This can potentially lead to 
lower losses and errors in some applications, i.e. some quantum information protocols. Our 
results could be also interesting in the field of nanophotonics where, in some applications one 
needs to produce light which closely resembles the symmetry of a certain system. It is likely 
that our method allows to introduce new degrees of freedom in this kind of applications, 
considering that our equation ([7]) is valid beyond the paraxial approximation and as stated 
above, in particular can be used with highly focused beams. 

In conclusion, our work closes an existing gap between the descriptions of angular mo- 
mentum in paraxial and nonparaxial electromagnetic fields. We have demonstrated that a 
paraxial beam with 0AM / and polarization p can be decomposed into multipolar modes 
with a fixed component of the angular momentum m = l+p and different components of the 
total angular momentum j. We have provided analytical rules to calculate the weight of the 
superpositions for several interesting cases, thus allowing a control over the superposition of 
the multipolar fields. These rules could be of utility in fields where the interaction of light 
and matter has to be controlled with precision. 
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